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Abstract 

We propose a novel relay augmentation strategy for extending the lifetime of a certain class of 
wireless sensor networks. In this class sensors are located at fixed and pre-determined positions 
and all communication takes place via multi-hop paths in a fixed routing tree rooted at the base 
station. It is assumed that no accumulation of data takes place along the communication paths and 
that there is no restriction on where additional relays may be located. Under these assumptions 
the optimal extension of network lifetime is modelled as the Euclidean fc-bottleneck Steiner tree 
problem. Only two approximation algorithms for this NP-hard problem exist in the literature: a 
minimum spanning tree heuristic (MSTH) with performance ratio 2, and a probabilistic 3-regular 
hypergraph heuristic (3RHH) with performance ratio VS + e. We present a new iterative heuristic 
that incorporates MSTH and show via simulation that our algorithm performs better than MSTH 
in extending lifetime, and outperforms 3RHH in terms of efficiency. 

1 Introduction 

Wireless sensor networks (WSNs) consist of small sensing devices that can be readily deployed in 
diverse environments to form distributed wireless networks for collecting information in a robust and 
autonomous manner. Although early research was mainly motivated by potential military uses, there 
are now many other important applications such as fault detection, environmental habitat monitoring, 
irrigation and terrain monitoring (see, e.g., pQ, [2]). An example of the latter is the proposed use 
of WSNs to provide an early warning system for bushfires in Australia. Deploying smart sensors in 
strategically selected areas can lead to early detection and an increased likelihood of success in fire 
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extinguishing efforts. Other applications include pollution control, climate control in large buildings, 
and medical applications using implantable devices. 

In this paper we assume that an initial WSN deployment phase has already taken place and that the 
locations of the sensors are known. In the relay augmentation phase a bounded number of relays are 
deterministically introduced to the sensing field, whereupon a fixed routing tree rooted at the base 
station and spanning all nodes is constructed. The power levels of the sensors and relays are then 
adjusted so that every node is able to transmit data to its parent in the routing tree. We assume that 
all sensors and relays transmit data at the same constant rate throughout the lifetime of the network. 

In most WSNs the bulk of a sensor's energy consumption is attributable to communication. In turn, 
the energy consumed during data transmission is proportional to some power (usually between 2 and 
4) of the communication distance. Therefore, under the above assumptions, the time till first node 
failure will be determined by the length of the longest edge in the routing tree, which we call the bottle- 
neck edge. There are many important WSN applications satisfying this property, including biomedical 
sensor networks [3], certain environmental monitoring networks, and natural disaster early detection 
systems. In many types of habitat monitoring applications, for example, the sensors are deployed at the 
locations of pre-determined data sources [2] . Sparsity of the sensor locations may then necessitate the 
introduction of relays to obtain a connected network and to optimize the lifetime of the resulting sparse 
network [4]. Furthermore, many WSN applications either utilise simple aggregation functions based on 
maximum/minimum or average readings, or use advanced data fusion techniques to limit accumulation 
[5], [6], [7J. Therefore our assumption of a uniform data transmission rate for all nodes is also realistic. 

The problem of optimal relay augmentation for extending the lifetime (i.e., time till first node-death) 
of WSNs of the above type has previously been modelled as the k-bottleneck Steiner tree problem (k- 
BST problem) or its dual, the minimum Steiner point tree problem [5].[5].[ID ] .[11 | . The objective of 
the fc-BST problem is, given a set of planar nodes (representing the sensors and called terminals), to 
introduce at most k new nodes (representing the relays and called Steiner points) so that the resulting 
minimum spanning tree on the complete set of nodes minimises the length of the bottleneck edge. As 
an example illustrating this problem consider Figure [T] where sensors are represented by black filled 
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circles, relays by open circles, and bottleneck edges by grey lines. In Figure (Ua) no relays are added 
and the given sensors are connected by a minimum spanning tree (MST). In Figure QJb) two relays are 
included, which clearly gives a solution with shorter bottleneck edge than in Figure [Ha). In fact, the 
solution in Figure D^b) is optimal for the given nodes when at most two Steiner points are allowed (note 
that other optimal solutions for the given terminals also exist). 

In this paper we introduce a new heuristic for this NP-hard problem, and provide simulation results 
that demonstrate a significant improvement in performance of our algorithm over the currently best 
performing algorithms. In particular, we compare our algorithm to the minimum spanning tree heuris- 
tic (MSTH) of Wang and Du [12]. Apart from making use of an exact algorithm for the 1-BST 
problem developed by Bae et al.[13|. our algorithm also incorporates MSTH as a module. In Section 
[5] we look at related work. Section [3] formalises our network model and describes some properties of 
optimal solutions to the fc-BST problem. In Section 0] we provide necessary details of existing heuristic 
algorithms for solving the fc-BST problem, as well as our new iterative 1-BST heuristic algorithm, and 
provide an analysis, with examples, of this algorithm in Section [5] Section [5] is devoted to an empirical 
demonstration of the performance of our algorithm. 

2 Related work 

There are several constraints and network performance metrics - for instance coverage, power usage, 
lifetime and cost - that predetermine an optimal WSN deployment strategy. The performance metric 
that we are interested in is lifetime since it is pertinent to such a large class of WSN scenarios. We 
define the lifetime of a network as the time till first node failure, which is generally due to battery 
depletion. 

Sensors are often overtaxed due to an uneven distribution of traffic flow. Many lifetime extending 
strategies in the literature therefore strive to dynamically adjust the routing topology in order to 
relieve the burden on these sensors; see [14j.|15j.|16j. Unfortunately, these techniques are only germane 
if multiple available paths exist between nodes in the network. In WSNs with low connectivity or where 
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routing protocols are limited by other application-specific or topological factors, deployment based 
strategies may be more valuable. 

Often the most highly burdened sensors are within close proximity to the base station [T7], hence 
the well-investigated objective of "doughnut" or "energy-hole" mitigation around the base station. 
Deploying additional relays close to the base station is one possible solution to the energy-hole problem, 
and this method has an analogue in density- varying random deployment strategies |18j . This approach, 
however, is only appropriate where there is a significant accumulation of data close to the base station. 
Given the assumption of a uniform data transmission rate for the networks dealt with in the current 
paper, these "energy-hole algorithms" are not comparable to our algorithm. 

Augmenting a network by deploying additional (non-sensing) relays is a powerful and relatively inex- 
pensive method of optimising many topology dependent WSN objectives. This concept is not new to the 
literature, and has been considered under various network models and objectives. Relay augmentation 
strategies exist that supplement the primary objective of lifetime extension with objectives related to 
coverage [19], connectivity j4] and balanced traffic flow [20] . [2T] . 

As mentioned above, here we model optimal lifetime extension as a fc-bottleneck Steiner problem. 
This problem was introduced by Saraffzedeh and Wong [22] . and various authors have considered it 
in the context of facility location science [23], [jH]. The problem is NP-hard; in fact, unless P=NP, no 
polynomial-time algorithm exists for the problem in the Euclidean plane with a performance ratio less 
than -\/2, where this ratio is defined as the theoretically largest possible value attained when dividing 
the length of the bottleneck edge produced by a given polynomial-time algorithm by the length of the 
bottleneck edge in an optimal solution. Wang and Du, the authors who demonstrated this complexity 
bound, present in 12] the first deterministic approximation algorithm for the fc-bottleneck problem in 
both the rectilinear and Euclidean planes. They also prove that the performance ratio of their algorithm, 
the afore-mentioned minimum spanning tree heuristic, is bounded above by 2 for each of these metrics. 
Du et al. |25] describe a probabilistic 3-regular hypergraph heuristic (3RHH) with performance ratio 
•\/3 + e which runs in - x poly(rt, k) time, where n is the number of given terminals. Recently Bae et 
al. [26], [13] developed exact algorithms with exponential complexity. 



4 



3 Sensor Network Model and Properties 



The energy consumption for each node in a WSN can be modelled, as in [21] , by 

E(r) = r a + c 

where r is packet transmission distance, c is a constant and is small relative to r Q , and a € [2,4]. The 
power consumption attributable to packet receipt is assumed to be negligible. Let l m nx(T) denote the 
length of the longest edge in a tree T. For any set X of n embedded terminals, let C(X, k) be the set 
of all trees interconnecting X and at most k other points in the plane. 

Definition. The k-bottleneck Steiner tree problem, or fc-BST problem, has the following form: 

Given A set X of n points (terminals) in M 2 , and a positive integer k. 

Find A set S (the Steiner points) of at most k points in M 2 , and a spanning tree TonlUS 
such that Cax(T) < 4iax(T') for any T' E C(X, k). 

We refer to T as a minimum k-bottleneck Steiner tree (k-MBST). In general, fc-MBSTs are not easy 
to construct. The following results, culminating in Corollary [3J allow us to restrict the search-space 
somewhat. The next lemma from [27] is an extension of the Swapping Algorithm found in [28) . It is a 
consequence of the matroid properties of minimum spanning trees. 

Lemma 1 Let T be a minimum spanning tree on the terminal set X , and let T be a spanning tree for 
X . We can transform T' to T by a series of edge swaps, where each swap involves replacing an edge 
ej G E(T') by ej € E(T) such that | e» | > \ej\, and at each stage the graph is a spanning tree. 

Corollary 2 If T' is a k-bottleneck Steiner minimum tree on X with Steiner points S then every 
minimum spanning tree on X U S is a k-bottleneck Steiner minimum tree on X . 

Proof. Let T be a minimum spanning tree on X U S. By Lemma [TJ we can transform T' to T by a 
series of edge swaps, each of which replaces an edge with another of the same length or shorter. ■ 



■5 



It is well-known in the literature that for any set of planar nodes an Euclidean MST exists with highest 
degree (i.e., the number of edges meeting at a node) at most 5. We therefore have the following result. 

Corollary 3 For any set of terminals there exists a k-MBST that is an MST on its complete set of 
nodes and such that the maximum node degree is at most 5. 

As recognized in [13) . the smallest spanning circle concept is fundamental to the fc-BST problem. A 
smallest spanning circle for a given set of nodes is a circle of smallest radius that includes all nodes 
on or within its boundary. The centre of a smallest spanning circle minimises the largest distance to 
a node of the given set. Therefore we may assume that every Steiner point in a fc-MBST is located at 
the centre of the smallest spanning disk of its neighbours (for if not then we could relocate the Steiner 
point to this location without increasing the length of the bottleneck); see Figure O 

The above properties now provide us with a simple (albeit time expensive) exact solution to the 1-BST 
problem: we select the cheapest tree after iterating through all subsets of at most five neighbours for the 
Steiner point and, for each iteration, locating the Steiner point at the centre of the smallest spanning 
disk on the subset and then constructing an MST on all nodes. The total complexity of this algorithm 
is therefore 0(n 6 logn). Fortunately the 1-BST can be solved in O(nlogn) time by using the algorithm 
of Bae et al. [13]. 

4 Algorithms 

Before providing the details of our main algorithm we briefly describe the algorithms MSTH and 3RHH. 
These algorithms will be compared to ours via simulation in Section [6] 

4.1 The Minimum Spanning Tree Heuristic 

The only deterministic algorithm in the literature for approximating a fc-MBST is the minimum spanning 
tree heuristic (MSTH) of Wang and Du [T^] . The authors prove that the performance ratio of MSTH 
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lies inclusively between \/3 and 2, and they conjecture that, if P ^ NP, y/3 is the best possible perfor- 
mance ratio of any polynomial-time algorithm for the Euclidean fc-BST problem. The time-complexity 
of MSTH is O(nlogn). In MSTH all Steiner points are of degree two and are referred to as beads. Let 
T be an MST on any q nodes. We define beading an edge e% of T with rij beads as placing equally 
spaced beads along a. Our main algorithm makes extensive use of MSTH as a "look-ahead" component, 

which we denote by BEAD. Subroutine BEAD is virtually identical to MSTH, except that in MSTH an 

| e .| 

MST is not part of the input (only the locations of the terminals are). For each e, in T let /(e,) = — . 

Ui + 1 

Algorithm 1 Subroutine BEAD (MSTH) 

Input: An MST T on q nodes and a positive integer j 

Output: A beaded MST T' 

1: Clear all beads of T, i.e., set n.; = for any i 

2: ex, e g _i be the edges of T 

3: Compute l{e,i) for each a 

4: Sort the edges in non-decreasing order of l(-) 

5: Add a bead to e; of largest l(-) value 

6: Update l(ei) 

7: Relocate the beads on a so that they are equally spaced 

8: Reset e,-'s position in the ordering 

9: Repeat Steps 5-8 until j beads have been added 

It is demonstrated in [T^] that Algorithm 1 is the best possible solution when all Steiner points are of 
degree two. 

4.2 The 3-Regular Hypergraph Heuristic 

In this section we discuss the only other heuristic for the fc-BST problem, namely the randomized 
algorithm of Du et al. |25j . The core of their algorithm is based on a randomized approximation 
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scheme for finding minimum spanning trees in weighted 3-uniform hypergraphs discovered by Promel 
and Steger [29]. This MST component of 3RHH, say MST3H, is the least efficient part of Du et al.'s 
algorithm, and we will need to discuss certain technical aspects of MST3H in order to demonstrate the 
superiority of our algorithm over 3RHH. 

As input to MST3H, a weighted 3-regular hypergraph containing 0(n 3 ) edges and with edge weights 
at most k is constructed by 3RHH. The edge-weights are then randomized by MST3H to a size of 
0{mnk), where m is the number of edges of the hypergraph. A matrix A is then constructed containing 
entries of order 0(2 W ), where w is some edge- weight. The determinant of A is calculated and is of size 
0((2 w ) n ). The determinant is therefore of order 2°^^. To find an MST in the given hypergraph the 
largest w such that 2 2w divides the determinant must be calculated. Due to the extremely large size of 
the determinant this takes very long in practice. 

As an example we implemented an instance with n = 10 and k = 3. Then 3-RHH produced a 
hypergraph with a determinant containing 496,631 digits, which made it practically intractable to find 
the relevant divisor as stipulated by the algorithm. The size of the determinant in this example is 
typical for instances with n = 10. 

4.3 A Naive Iterative 1-Bottleneck Steiner Tree Heuristic 

A naive version of the iterative 1-BST heuristic (or Il-BSTH) on a set X of terminals might proceed 
as follows. An optimal 1-BST on X is constructed by Bae et al.'s algorithm. This sequence of routines 
is repeated k — 1 times, i.e., until k new points have been added. 

In order to compare the performance of this heuristic against MSTH we implemented both in C (with 
minor modifications as described below in Subroutine CONVERT), and tested them on uniformly 
distributed terminal sets. Roughly speaking, simulations show a maximum average improvement of 
approximately 4% in maximum edge length for naive Il-BSTH over MSTH. This average can be in- 
creased to about 5% by combining naive Il-BSTH and MSTH into a meta-heuristic, which shows that 
MSTH can outperform naive Il-BSTH on some terminal sets. A simple example where this occurs is 



described next, and is illustrated in Figures [3] and 21 in these figures (and throughout) the solid circles 
are terminals and the open circles are Steiner points. 

Example 1 Letk = 2 and consider four terminals with coordinates ti = (2.0,9.1), t% = (3.0,8.6), t3 = 
(4.6,3.1), t 4 = (8.6,9.2). Then pi = 2.85 and P2 = 3.63 represent the longest edges in the solutions of 
MSTH and naive Il-BSTH respectively. 

Fig. [4] also illustrates one of the modifications made to naive Il-BSTH, namely that of correcting any 
bent edges that may appear in the final solution; a bent edge actually consists of two edges incident, 
at an angle of less than 180°, to a Steiner point of degree two. 



4.4 The Pre-beaded Iterative 1-Bottleneck Steiner Tree Heuristic 



Fig. [4] highlights a major draw-back of iterative greedy heuristics: an optimal solution at some itera- 
tion may undercut the gains of subsequent iterations. We notice that, bent-edge correction aside, the 



tree from Fig. 4(c) provides no decrease in maximum edge-length over Fig. 4(a) This shortcoming 



permeates the class of solutions constructed by naive Il-BSTH. A look-ahead routine is therefore war- 
ranted as a modification to Il-BSTH, but only if gains in performance are made by such a routine over 
the afore-mentioned meta- heuristic. The pre-beaded iterative 1-BST heuristic (pre-beaded Il-BSTH), 
described in Algorithm 2, utilises subroutine BEAD for this purpose. 



Algorithm 2 Pre-beaded Il-BSTH 

Input: A set X of n points in the plane, and a positive integer k 

Output: An approximate fc-MBST T interconnecting X and k new points 

1: Construct a minimum spanning tree To on X 

2: Let p = c = 

3: while c < k and p < k 

4: Run BEAD with input T c and k — 1 — p, and output T c ' 

5: Run Bae's algorithm for 1-BST on T c ' with output T' and s 
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6: Run CONVERT with input T' and s, and output T" 

7: Let p' be the number of non-bead Steiner points of T' 

8: Run BEAD with input T" and k — p', and output T* 

9: end for 

10: Increment c by 1, let T c be the cheapest T* produced 

11: Let p be the number of non-bead Steiner points of T c 

12: end while 



Pre-beadcd Il-BSTH also contains the procedure CONVERT, which consists exactly of the modifica- 
tions that were made to naive Il-BSTH for its implementation in C. 

Algorithm 3 Subroutine CONVERT 

Input: A tree T interconnecting a set of terminals X and a set of Steiner points S, and an element 
s g 5" 

Output: The modified tree T' 

1: Delete all degree-one Steiner points of T 

2: Correct all bent edges and relocate beads to be equally spaced along edges 

3: Relocate s at the Chebyshev point of its neighbours, and add the edges incident to s 



5 Analysis of Pre-beaded Il-BSTH 

There are three aspects of pre-beaded Il-BSTH that we examine more closely in this section. The first 
is the look-ahead component actualised by the call to BEAD in Line 12 of Algorithm 2. Calculating 
an optimal sequence of Steiner point additions can generally not be done greedily, and calculating 
every potential sequence would take us into exponential complexity. However, we can remain within 
polynomial complexity if we assume, at any given iteration of the algorithm, that the remaining sequence 
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of additions will consist of beads only. This assumption can be viewed as providing a lower bound for 
the overall (including past and future iterations) performance of an optimal Steiner point addition 
sequence. This lower bound is clearly equal to the performance of MSTH, and hence we are guaranteed 
that pre-beaded Il-BSTH never performs worse than MSTH. This result on its own is not satisfactory, 
since the meta heuristic combining the naive version of Il-BSTH and MSTH has the same property. 
But as will be demonstrated in the next section, algorithm pre-beaded Il-BSTH performs significantly 
better than the meta heuristic. 

We use Fig. [5] to illustrate the benefit of the look-ahead component: we compare the solutions 
constructed by naive Il-BSTH and pre-beaded Il-BSTH on a set of six terminals ti = (968.4,506.4), 
t 2 = (3.9, 86.8), t 3 = (188.8, 7.5), t 4 = (779.2, 675.9), t 5 = (238.1, 644.4), t 6 = (620.6, 2.4), and where 
k = 2. The three longest edges of the MST on these terminals are (in descending order of length) 



ei,e 2 ,e 3 ; see Fig. 5(a 



The naive algorithm, Fig. 5(b) selects {ti, t 4 , ts, te} as the optimal neighbour set for the first Steiner 
point, si. The Chebyshev centre of this set is defined by the subset {ti,t5,te}, and the subsequent 
tree, after updating the MST to include s±, has e?, as its longest edge. The second iteration's optimal 
point happens to be a bead subdividing e^, and the longest edge of the final tree is therefore = s±te 



(equivalently Sits or Siti) and has a length of 389.87 units. Fig. 5(c) shows the final tree produced by 
pre-beaded Il-BSTH. Here one may view the algorithm as initially adding a single bead to each of edges 
ei and e 2 ; it then optimally repositions the bead on ei to create the Steiner point with neighbour-set 
{t3,t5,t6j-. Finally it looks for an improved position for the bead on e 2 , but determines that that 
bead is already an optimal Steiner point with respect to the other nodes. The resulting longest edge is 
e§, with | eg ] = 374.46. Note also that naive Il-BSTH greedily deletes the two longest edges, ei,e2, in 
its first iteration, resulting in a final tree with relatively long edges incident to the first Steiner point. 
Pre-beaded Il-BSTH, however, deletes edges e\ and e$ in the first iteration because it assumes that e 2 
will at least be beaded in the second iteration. The result is a Steiner point with relatively shorter 
incident edges. This example is simple, but typical of many cases where the pre-beaded algorithm 
outperforms the naive version. 
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The second aspect we would like to examine concerns the call to BEAD in Line 4. This call ensures 
that there are k — 1 Steiner points in the input tree before any main iteration. Consequently our 
algorithm may be viewed as beading the entire initial MST, and then one-by-one removing beads and 
adding them back optimally. It is possible to construct a version of Il-BSTH that does not utilise this 
BEAD call but still contains the look-ahead routine. This so called Post-beaded Il-BSTH performs very 
similarly to our pre-beaded algorithm, but falls short in certain instances, of which Fig. [5] is a case in 
point. 



Before the first iteration of pre-beaded Il-BSTH the MST in Fig. 6(a) already contains a bead on the 
longest edge tit4. As opposed to the post-beaded version, our algorithm is able to utilise this bead as 
a potential neighbour for the Steiner point Si, and the result is the simultaneous introduction of two 
Steiner points during the first iteration. 



6 Empirical Results 

Multiple /c-BST problem instances, uniformly distributed (i.e., randomly with a uniform distribution) 
in a 10,000 by 10,000 point region, were generated and solved in MATLAB for MSTH, post-beaded 
Il-BSTH, and pre-beaded Il-BSTH. We focussed on terminal sets of size n = 600 and k = 20, 40, 200 
and performed 300 simulations for each case. Tables HHSl and Figures [7] and [9] present the simulation 
results. We also performed 300 simulations for instances with terminals distributed non-uniformly in 
a 10,000 by 10,000 point region, with results depicted in Figured] (once again we focussed on n = 600 
and k = 20, 40, 200). In particular, for this last group of simulations we selected a random point to 
represent a base station and then distributed the points with an increasing density towards the base 
station. The run times for our algorithm are not as important as algorithms which are dynamic in 
nature (such as for instance routing protocols). However, it is not unreasonable to envision sensor 
network applications arising in the future that have many more terminals than just a few thousand; 
in other words we must still show that our algorithms scale well. We have therefore also included run 
time results in our simulations. In Table [4] we present average run time results for 3RHH on instances 
of up to 6 terminals and 3 Steiner points. 
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Table 1: MSTH 



k 


Ave. lifetime (s) a = 2 


Ave. lifetime (s) a = 4 


Ave. run time (s) 


20 


19.25 


0.037 


1.05 


40 


20.12 


0.040 


1.12 


60 


21.37 


0.046 


1.28 


80 


23.74 


0.056 


1.29 


100 


24.65 


0.061 


1.33 


120 


25.68 


0.066 


1.44 


140 


26.64 


0.071 


1.47 


160 


28.50 


0.081 


1.51 


180 


30.95 


0.096 


1.58 


200 


31.89 


0.102 


1.68 



Tabic 2: Post-beaded Il-BSTH 



k 


Avg. lifetime (s) a = 2 


Avg. lifetime (s) a = 4 


Avg. run time (s) 


20 


20.53 


0.042 


1.88 


40 


22.19 


0.049 


2.31 


60 


24.31 


0.059 


2.97 


80 


27.11 


0.073 


3.25 


100 


28.58 


0.082 


3.88 


120 


30.44 


0.092 


4.41 


140 


31.06 


0.096 


4.89 


160 


35.66 


0.127 


5.39 


180 


38.63 


0.149 


5.81 


200 


40.55 


0.164 


6.33 
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Table 3: Pre-beaded Il-BSTH 



k 


Avg. lifetime (s) a = 2 


Avg. lifetime (s) a = 4 


Avg. run time (s) 


20 


21.81 


0.048 


1.98 


40 


21.56 


0.047 


2.22 


60 


26.34 


0.069 


3.02 


80 


28.20 


0.080 


3.20 


100 


29.08 


0.085 


3.89 


120 


30.88 


0.095 


4.49 


140 


32.77 


0.107 


4.77 


160 


35.37 


0.125 


5.29 


180 


40.24 


0.162 


6.12 


200 


44.35 


0.197 


6.56 



Table 4: Run times for 3RHH 



n 


k 


Avg. run time (s) 


3 


2 


23.52 


4 


2 


180.41 


5 


3 


7653.23 


6 


3 


28845.49 



In our simulations we observed that the performance of pre-beaded Il-BSTH relative to MSTH roughly 
increases from k = 0, reaching a peak in the approximate range 0.2 < ^ < 0.34, whereafter performance 
decreases at a decreasing rate. We also observed near-asymptotic behaviour as k increases beyond k = n 
(see Figure [T0l where ^bsth an< ^ ^msth are defined below). The limiting behaviour might not continue 
indefinitely, but we would like to state the following lemma and corollary. 

For any set X of n terminals let £q PT (X), ^msth(^)j an d ^bsthO^O represent the length of the 
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longest edge in three potential solutions to the fc-BST problem on X; namely the optimal solution, 
the solution generated by MSTH, and the solution generated by pre-beaded Il-BSTH respectively. Let 
Lsmt(X) be the total length of the classical Steiner minimal tree (SMT) on X (i.e., the shortest total 
length tree spanning X where any number of Steiner points may be introduced), and Lmst(X) be 
the total length of the MST on X. We let £(X, typel, type2) = Jim (£% ypel (X) / 1 £ k typc2 {X)) , where 
typel,type2 G {OPT, MSTH, BSTH}. 

Lemma 4 £(X, MSTH, OPT) = L MS t(X)/ L SM t(X). 

Proof. As k — >• oo, the length of every edge in the optimal beading of the MST on X approaches 
^msth(^)' an d hence ^msthG^O — ^ -^mst(A). Similarly, as k — > co in the optimal bottleneck tree, at 
most n — 2 Steiner points have degree greater than 2, and hence almost all Steiner points are beads. So 
again k.£Q PT (X) approaches the total length of the minimum length tree interconnecting X , namely 
the SMT on X. ■ 

Let E(-) be the expected value of some parameter determined by X, where X is uniformly distributed, 
and let Eo = E (XmstPO / Lsmt(X)) ■ Gilbert and Pollack [3D] provide an upper-bound of 1.064 for Eo. 

Corollary 5 If E (£(X, BSTH, MSTH)) =p then 
E (l(X, BSTH, OPT)) = p.E a 

Proof. 

E(^(X,BSTH, OPT)) 

= E (£(X, MSTH, OPT) • £(X, BSTH, MSTH)) 
= E (£(X, MSTH, OPT)) • E (£(X, BSTH, MSTH)) 
= E(L MST (X)/L SMT (X)) ■ p. m 

Figure HOl suggests that E (£(X, BSTH, MSTH)) is approximately 0.985. Based on this value, and by 
using the above corollary, we can estimate the performance ratio of Il-BSTH on uniformly distributed 
sets as k tends to infinity: E (£(X, BSTH, OPT)) = p.E 0.985 x E < 1.04804, where the final 
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inequality comes from the above mentioned upper bound on Eo given in |30| . 

7 Conclusion 

One can model the problem of the optimal lifetime extension of homogeneous WSNs as the Euclidean 
/c-bottleneck Steiner tree problem. In this paper we propose a new approximation algorithm for the 
fc-BST problem, namely the pre-beaded iterative 1-BST heuristic, and test its performance relative to 
the minimum spanning tree heuristic (MSTH) of Wang and Du, which is currently the best performing 
heuristic for the problem. By incorporating MSTH as a "look-ahead" component and utilising Bae 
et al.'s algorithm for solving the 1-bottleneck Steiner tree problem, we develop an algorithm that 
significantly outperforms MSTH for uniformly distributed terminal sets. We also demonstrate that the 
probabilistic 3-regular hypergraph heuristic (3RHH) of Du et al. is impractical even for small n. To 
our knowledge this is the first well-performing tractable approach to solving the fc-BST problem. 
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(a) A minimum spanning tree on six 
given nodes (sensors) 



(b) A spanning tree on the given nodes 
which includes two Steiner points (re- 
lays) 



Figure 1: Two potential solutions to an instance of the 2-bottleneck Steiner tree problem 
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(a) 



(b) 



Figure 2: Steiner points may be relocated, without increasing the bottleneck, so that they lie at the 
centre of the smallest spanning circle of their neighbours 
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(a) Four given terminals (b) Longest edges are 

plus their MST beaded 

Figure 3: The solution constructed by MSTH as described in Example Q] 
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(a) First Steiner point (b) Second Steiner (c) Tree is updated (d) Bent edge is cor- 

is added point is added rected 

Figure 4: The solution constructed by naive Il-BSTH as described in Example Q] 
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(a) MST 



(b) Naive Il-BSTH solu- (c) The pre-beaded II- 

tion BSTH solution 



Figure 5: The benefit of looking ahead 
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(a) MST 



(b) Pre-beaded Il-BSTH 
solution 



Figure 6: The benefit of pre-beading 
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Figure 7: Comparison of the lifetime of networks constructed by Pre-beaded Il-BSTH, Post-beaded 
Il-BSTH, and MSTH for a uniform distribution when a — 4 
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Figure 8: Comparison of the lifetime of networks constructed by Pre-beaded Il-BSTH, Post-beaded 
Il-BSTH, and MSTH for a non-uniform distribution when a = 4 
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Figure 9: Comparison of the running times for network construction in Pre-beaded Il-BSTH, Post- 
beaded Il-BSTH, and MSTH 
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Figure 10: The ratio 4sth/^msth for n = 200 and 400 < k < 1000 
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